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Hie  wide  availability  of  electronic  computers  today  allows 
calculations  of  a  detail  and  accuracy  which  was  impossible  a 
few  years  ago,  but  these  computer  calculations  are  only  as  valid 
as  the  input  data  upon  which  they  are  based.  In  the  case  of 
ships’  hull  structures,  the  increased  calculational  capability 
has  meant  that  loads  acting  on  the  hull  must  be  known  more  accu¬ 
rately  than  ever  before. 

A  major  portion  of  the  effort  of  the  Ship  Structure  Committee 
research  program  has  been  devoted  to  improving  capability  of 
determining  hull  loads.  This  report  and  the  two  which  follow  it 
concern  a  project  directed  towards  this  end,  which  involved  the 
development  of  a  computer  program  to  calculate  these  loads. 

This  report  contains  a  description  of  the  development  and 
verification  of  the  program  for  predicting  hull  loads.  SSC-230, 
Program  SCORES- -Ship  Structural  Response  in  Waves,  contains  the 
details  of  the  computer  program  and  SSC-231,  Further  Studies  of 
Computer  Simulation  of  Slamming  and  Other  Wave- Induced  Vibratory 
Structural  Loadings  on  Ships  in  Waves,  contains  further  details 
on  the  use  of  the  analysis  method  for  prediction  of  other  hull 
loadings . 

Comments  on  this  report  would  be  welcomed. 
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An  analytical  method  for  the  determination  of  conventional  merchant 
ship  motions  and  wave-induced  moments  in  a  seaway  is  developed.  Both  verti¬ 
cal  and  lateral  plane  motions  and  loads  are  considered  for  a  ship  travelling 
at  any  heading  in  regular  waves  and  in  irregular  long  or  short  crested  seas. 
Strip  theory  is  used  and  each  ship  hull  cross-section  is  assumed  to  be  of 
Lewis  form  shape  for  the  purpose  of  calculating  hydrodynamic  added  mass  and 
damping  forces  in  vertical,  lateral  and  rolling  oscillation  modes.  The 
coupled  equations  of  motion  are  linear,  and  the  superposition  principle  is 
used  for  statistical  response  calculations  in  irregular  seas.  All  three  pri¬ 
mary  ship  hull  loadings  are  determined,  i.e.  vertical  bending,  lateral  bend¬ 
ing  and  torsional  moments,  as  well  as  shear  forces,  at  any  point  along  the 
length,  with  these  responses  only  representing  the  low  frequency  slowly  vary¬ 
ing  wave  loads  directly  induced  by  the  waves. 


A  computer  program  that  carries  out  the  calculations  was  developed, 
and  is  fully  documented  separately.  The  results  of  the  method  are  evaluated 
by  comparison  with  a  large  body  of  model  test  data.  The  comparison  extends 
over  a  wide  range  of  ship  speeds,  wave  angles,  wave  lengths,  and  loading 
conditions,  as  well  as  hull  forms.  The  agreement  between  the  calculations 
and  experimental  data  is  generally  very  good.  Thus,  a  method  is  available 
for  use  in  the  rational  design  of  the  ship  hull  main  girder  structure. 
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ABSTRACT 


An  analytical  method  for  the  determination  of  conventional  merchant 
ship  motions  and  wave-induced  moments  in  a  seaway  is  developed.  Both  verti¬ 
cal  and  lateral  plane  motions  and  loads  are  considered  for  a  ship  travelling 
at  any  heading  in  regular  waves  and  in  irregular  long  or  short  crested  seas. 
Strip  theory  is  used  and  each  ship  hull  cross-section  is  assumed  to  be  of 
Lewis  form  shape  for  the  purpose  of  calculating  hydrodynamic  added  mass  and 
damping  forces  in  vertical,  lateral  and  rolling  oscillation  modes.  The 
coupled  equations  of  motion  are  linear,  and  the  superposition  principle  is 
used  for  statistical  response  calculations  in  irregular  seas.  All  three  pri¬ 
mary  ship  hull  loadings  are  determined,  i.e.  vertical  bending,  lateral  bend¬ 
ing  and  torsional  moments,  as  well  as  shear  forces,  at  any  point  along  the 
length,  with  these  responses  only  representing  the  low  frequency  slowly  vary¬ 
ing  wave  loads  directly  induced  by  the  waves. 

A  computer  program  that  carries  out  the  calculations  was  developed, 
and  is  fully  documented  separately.  The  results  of  the  method  are  evaluated 
by  comparison  with  a  large  body  of  model  test  data.  The  comparison  extends 
over  a  wide  range  of  ship  speeds,  wave  angles,  wave  lengths,  and  loading 
conditions,  as  well  as  hull  forms.  The  agreement  between  the  calculations 
and  experimental  data  is  generally  very  good.  Thus,  a  method  is  available 
for  use  in  the  rational  design  of  the  ship  hull  main  girder  structure. 
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NOMENCLATURE 


a  wave  amplitude 

a' ,b,c* ,d,e,g'  =  coefficients  in  vertical  (heave)  equation  of 

motion 

a^_j  =  ;ceff icients  in  lateral  plane  equations  of  motion 
a”7  =  mean  squared  response  amplitude 

A,B,C,D,E,G '  =  coefficients  in  vertical  plane  (pitch)  equation 

of  motion 

A  =  ratio  of  generated  wave  to  heave  amplitude  for  vertical 

motion-induced  waves 

A^3  =  sectional  vertical  added  mass 

A,  B  =  coefficients  in  two-parameter  spectrum  equation 
B*  =  local  waterline  beam 

B*  =  waterline  beam  amidships 


=  lateral  bending  moment 

=  vertical  bending  moment 

=  wave  speed  (celerity) 

=  local  section  area  coefficient 

=  total  local  lateral  loading  on  ship 

=  total  local  vertical  loading  on  ship 

=  total  local  torsional  loading  on  ship 

=  sectional  hydrodynamic  moment,  about  x  axis,  on  ship 

=  sectional  lateral  hydrodynamic  force  on  ship 

=  sectional  vertical  hydrodynamic  and  hydrostatic  force 
on  ship 

=  Froude  number 

=  sectional  lateral  added  mass  due  to  roll  motion 

=  acceleration  of  gravity 
=  center  of  gravity  of  ship 
=  initial  metacentric  height  of  ship 
=  mean  section  draft 
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H  =  sectional  draft 

Hjyj  =  significant  wave  height 

I  =  sectional  added  mass  moment  of  inertia 

1,1,1  =  mass  moments  of  inertia  of  ship  about  x,y,z  axes 

^  respectively 

I  =  mass  product  of  inertia  of  ship  in  x-z  plane 

X  z 

k  =  wave  number 

Kw  =  wave  excitation  moment,  about  x  axis,  on  ship 

L  =  ship  length 

m  =  mass  of  ship 

M  =  sectional  lateral  added  mass 

s 

M  =  wave  excitation  moment,  about  y  axis,  on  ship 

M  =  sectional  added  mass  moment  of  inertia  due  to 

s*  lateral  motion 

A  =  sectional  roll  damping  moment  coefficient  due  to 

wave  effects 

N*  =  sectional  roll  damping  moment  coefficient  due  to 
r  viscous  and  bilge  keel  effects 

Ng  =  sectional  lateral  damping  force  coefficient 

Nw  =  wave  excitatio..  moment,  about  z  axis,  on  ship 

N'  =  sectional  vertical  damping  force  coefficient 

N  =  sectional  lateral  damping  force  coefficient  due  to 

roll  motion 

N  ,  =  sectional  damping  moment  coefficient  due  to  lateral 

s*  motion 

OG  =  vertical  distance  between  waterline  and  center  of 
gravity,  positive  up 

S  =  local  section  area 

S (ui,y)  =  directional  spectrum  of  the  seaway 

S^(u,y)=  response  spectrum,  for  a  particular  response 

=  frequency  spectrum 
S  2 ( u )  =  spreading  function 

t  =  time 

$  =  mean  wave  period 

T,,T2,...  =  coefficients  in  lateral  plane  wave  excitation 

equations 

T^(w,y)=  response  amplitude  operator 
TMX  =  torsional  moment 

U  =  wind  speed 

vw  -  lateral  orbital  wave  velocity 

V  =  ship  forward  speed 
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=  horizontal  axis  in  direction  of  forward  motion  of 
ship  (along  length  of  ship) 

=  axis  fixed  in  space 

=  location  along  ship  length  at  which  moments  are  determined 
=  x  coordinates  at  stern  and  bow  ends  of  ship,  respectively 
-  horizontal  axis  directed  to  starboard;  sway 
=  lateral  wave  excitation  force  on  ship 
=  vertical  axis  directed  downwards;  heave 

=  vertical  space  coordinate,  from  undisturbed  water  surface, 
positive  downwards 

=  sectional  center  of  buoyancy,  from  waterline 

=  vertical  wave  excitation  force  on  ship 

=  angle  between  wave  propagation  direction  and  ship 
forward  motion 

=  local  mass  gyradius  in  roll  (about  x  axis) 

K,a,^,T,u  =  phase  angles  (leads)  of  heave,  pitch,  vertical 

bending  moment,  sway,  yaw,  roll,  lateral  bending 
moment,  torsional  moment,  respectively 

=  local  mass 

=  local  vertical  center  of  gravity,  from  CG,  positive  down 
=  fraction  of  critical  roll  damping 

=  surface  wave  elevation,  positive  upwards  from  undisturbed 
water  surface 

=  pitch  angle,  positive  bow-up 
=  wave  length 

=  wave  direction  relative  to  predominant  direction 
=  density  of  water 

=  roll  angle,  positive  starboard-down 
=  yaw  angle,  positive  bow-starboard 
=  circular  wave  frequency 
=  circular  frequency  of  wave  encounter 
=  natural  roll  frequency 


Subscripts 


average  (statistical) 

amplitude 

root-mean-squared 

significant  (average  of  1/3  highest) 
average  of  1/10  highest 
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A 

INTRODUCTION 


In  order  to  investigate  the  utility  of  a  computer  simu¬ 
lation  approach  for  determining  ship  bending  moment  responses  in 
waves,  a  research  program  was  instituted  under  the  sponsorship 
of  the  Ship  Structure  Committee  with  the  aid  of  an  advisory  panel 
appointed  by  the  National  Academy  of  Sciences.  The  original 
program  was  considered  to  be  made  up  of  three  separate  phases  of 
work  which  include: 

1.  An  assembly  of  a  system  of  equations  that  would  adequately 
describe  ship  structural  responses  due  to  the  effects  of 
waves . 

2.  The  conversion  of  these  equations  to  a  computer  program  or 
to  the  design  of  a  computer  analog. 

3.  Computer  evaluation  of  the  ship  response  mathematical  model, 
with  the  verification  of  the  entire  procedure  provided  by 
such  an  evaluation. 

The  first  phase  of  this  work,  which  was  the  development 
of  a  mathematical  model,  was  completed  and  described  in  a  final 
report  [1] .  A  mathematical  model  was  developed  under  that  program, 
where  equations  for  determining  wave -induced  bending  moments  in 
the  vertical  and  lateral  planes  were  established.  In  addition,  a 
method  of  treatment  for  including  effects  due  to  slamming  was 
outlined,  where  the  occurrence  of  slamming  was  evidenced  by  "whip¬ 
ping"  responses  that  may  be  ascribed  to  nonlinear  forces  generated 
due  to  bow  flare. 

The  second  phase  of  work  in  this  program,  which  has  been 
completed  and  described  in  report  form  [2] ,  is  devoted  to  the 
conversion  of  the  equations  developed  in  Phase  I  into  a  computer 
program.  The  linearized  vertical  plane  motions  and  vertical  bend¬ 
ing  moment  response  operators  for  a  ship  were  determined  by  a 
digital  computer  program  for  the  case  of  head  seas,  and  this  program 
was  then  generalized  to  the  case  of  oblique  headings  between  the 
ship  and  seaway.  Modifications  of  the  basic  head  sea  program  have 
been  carried  out  under  Phase  II  of  this  overall  program  (see  [2]), 
and  hence  these  quantities  are  amenable  to  computation  by  a  digital 
computer.  Further  modifications  to  incorporate  a  given  wave  spectrum, 
together  with  a  directional  spreading  factor  to  account  for  shcrt- 
crestedness,  will  allow  this  program  to  compute  the  power  spectra 
of  vertical  bending  moments  on  a  ship  in  irregular  short-crested 
seas. 

Since  lateral  bending  moments  occur  in  oblique  sea  conditions, 
and  since  they  have  significant  magnitude  in  certain  cases  relative 
to  the  vertical  bending  moment  for  that  same  heading,  a  program  for 
computation  of  lateral  bending  moments  has  also  been  developed  as 
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well  in  Phase  II.  While  the  hydrodynamic  data  for  this  particular 
structural  component  is  not  extensively  treated  in  the  available 
literature  as  is  the  case  of  vertical  motions  and  structural 
responses,  there  is  sufficient  basic  information  that  allows  a 
similar  treatment  to  be  applied  to  the  lateral  loads  although 
no  computer  programs  to  calculate  the  sectional  added  mass  and 
damping  due  to  lateral  and  rolling  motions  had  been  established 
previously.  Thus  lateral  bending  moment  spectra  can  then  also  be 
obtained  for  a  particular  input  wave  spectrum,  and  these  results 
can  be  combined  with  those  for  the  vertical  bending  moment,  if 
desired. 

The  work  described  in  the  present  report  treats  the 
analytical  determination  of  one  aspect  of  sea  loads,  viz.  the 
determination  of  wave-induced  moments  that  are  slowly  varying  in 
time  and  have  the  same  frequency  characteristics  as  the  encountered 
waves.  Other  sea-induced  loadings,  such  as  whipping,  slamming 
and  springing,  which  are  of  higher  frequency,  must  not  be  neglected 
in  an  overall  design,  and  analytical  work  to  cover  tnese  subjects 
has  also  been  carried  out  under  the  present  contract,  which  will 
be  reported  separately.  Wave-induced  moments  depend  both  on  the 
motion  responses  of  the  ship  and  the  wave-excitation  loads  them¬ 
selves.  These  factors,  in  turn,  depend  on  the  ship  geometry  and 
mass  distribution,  as  well  as  on  the  particular  wave  conditions. 

The  present  report  is  a  continuation  of  work  previously 
reported  ir.  [1]  and  [2].  While  much  of  the  previous  analytical 
results  with  respect  to  wave-induced  moments  are  repeated  here, 
analysis  procedures  for  wave-induced  moments  are  slightly  expanded 
and  refined  in  the  present  report  as  well  as  extended  to  include 
torsional  moments,  and  the  results  of  more  extensive  computer 
calculations  based  on  these  procedures  are  evaluated  by  comparison 
with  experimental  data.  The  digital  computer  program  (SCORES) 
developed  in  the  course  of  this  work  is  fully  documented  separately 
13]. 

The  present  results  apply  to  conventional  merchant  ship 
hull  forms.  Consideration  is  given  in  the  analysis  to  both  vertical 
and  lateral  plane  motxon  responses  and  wave-induced  moments,  with 
the  ship  advancing  at  any  heading  with  respect  to  the  waves.  The 
wave  environment  can  be  represented  as  either  regular  sinusoidal 
waves,  a  long-crested  (unidirectional)  seaway  of  specific  spectral 
form  or  a  fully  short-crested  seaway,  using  various  wave  energy 
spectral  formulations.  The  three  primary  ship  hull  loadings  that 
are  considered  are,  vertical  bending  moment,  lateral  bending  moment, 
and  torsional  moment,  with  primary  emphasis  upon  vertical  and  lateral 
beiding  (the  related  shear  forces  are  also  determined  in  this  work) . 

Since  the  necessary  inputs  to  the  wave-induced  moment  deter¬ 
mination  are  the  rigid  body  ship  motion  responses,  these  must  be 
obtained  initially.  The  equations  of  ship  motion  are  taken  to  be 
linear  and  coupled  only  within  each  plane.  That  is,  heave  and 
pitch  motions  are  coupled  in  the  vertical  plane,  and  sway,  yaw, 
and  roll  motions  are  considered  coupled  in  the  lateral  plane.  The 
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equations  are  solved,  or  more  precisely  the  terms  in  the  equations 
are  computed  by  application  of  "strip"  theory,  where  local  forces 
on  each  ship  section,  or  strip,  are  evaluated  independently,  with¬ 
out  allowing  for  influence,  or  interaction,  among  sections.  This 
method  was  originally  derived  by  Korvin-Kroukovsky  [4] ,  and  in 
collaboration  with  Jacobs  [5] ,  for  vertical  plane  motions,  and  has 
subsequently  been  adopted  and  expanded  by  many  investigators. 

The  hydrodynamic  forces  at  each  station  which  enter  into  the 
equations  of  motion  are  obtained  by  a  potential  flow  solution  for 
an  equivalent  "Lewis"  form  section  shape  [6] .  In  general,  the 
Lewis  form  shape,  defined  simply  by  two  paramters  (beam-draft  ratio 
and  section  area  coefficient) ,  is  considered  to  be  a  fairly  close 
representation  of  section  shapes  found  in  conventional  merchant 
ship  hulls,  without  a  large  bulb  at  the  bow.  The  hydrpdynamic  forces, 
added  mass  and  damping,  are  obtained  for  vertical  section  oscillations 
by  the  method  developed  by  Grim  [7] ,  and  for  lateral  and  rolling 
oscillations  by  the  method  of  Tasai  [8]. 

The  present  work  is  aimed  at  verifying  the  capability  of  a 
digital  computer  technique  in  providing  valid  information  for 
evaluating  wave-induced  ship  structural  responses  under  various 
environmental  conditions,  for  ships  having  conventional  hull  forms. 
This  is  achieved  by  applying  the  method  of  computation  to  a  number 
of  particular  cases,  which  represent  computer  experiments  that 
point  out  simplications,  improvements,  etc.  that  can  be  incorporated 
in  a  final  computer  program.  The  program  will  provide  codification 
of  various  elemental  steps,  specific  subroutines  for  computing 
separate  items  such  as  sectional  hydrodynamic  forces,  etc.,  and  the 
computational  experiments  are  used  to  establish  a  final  formulation 
of  a  complete  and  efficient  digital  computer  program  that  will 
produce  structural  response  information  with  a  minimized  computer 
time  and  cost.  A  fully  documented  computer  program,  including  a 
description  of  data  input,  output  forms,  flow  charts,  and  the  pro¬ 
gram  listing  are  given  in  [3] .  The  results  of  extensive  computa¬ 
tions  for  a  number  of  ships,  for  which  model  test  data  are  avail¬ 
able,  are  presented  in  the  present  report  together  with. a  comparison 
between  the  computations  and  the  experiments. 


ANALYTICAL  METHOD 


The  basic  analytical  procedures  for  the  determination'  of  the 
wave-induced  moments  were  presented  originally  in  [1] .  In  the 
course  of  the  work,  certain  additions  and  modifications  to  the  . 
original  development  have  been  deemed  advisable.  Therefore,  the 
full  analytical  treatment  is  presented  below,  with  the  refinements 
included. 

The  coordinate  system  relationship  between  the  water  wave 
system  and  the  ship  coordinate  axes  is  shown  in  Figure  1.  Whereas 
in  the  previous  work,  separate  axes  conventions  were  employed  for 
the  vertical  and  lateral  motions  cases,  a  single  ship  axes 
coordinate  system  is  now  used.  All  the  equations  of  motions  are 
formulated  relative  to  a  right-handed  cartesian  coordinate  axes 
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system  whose  origin  is  located  at  the  center  of  gravity  of  the  ship, 
G,  and  with  the  x-axes  positive  toward  the  bow  (in  the  direction  of 
forward  motion),  the  y-axis  positive  to  starboard,  and  the  z-axis 
positive  downward.  These  axis  are  defined  to  have  a  fixed 
orientation,  i.e.  they  do  not  rotate  with  the  ship,  but  they  can 
translate  with  tiie  ship.  The  ship  angular  motions  are  considered 
to  be  small  oscillations  about  the  mean  position  defined  by  the 
axes. 

The  wave  propagation,  at  speed  c,  is  considered  fixed  in 
space.  The  ship  then  travels,  at  speed  V,  at  some  angle  8  with 
respect  to  the  wave  direction.  The  wave  velocity  potential,  for 
simple  deep-water  waves,  is  then  defined  by: 

4>  =  ace~^z  cos  k  (x'  +  ct)  (1) 

w 

where  a  =  wave  amplitude 
c  =  wave  speed 

,  ,  2tt 

k  =  wave  number  =  r—  =  — 

X  g 

X  =  wave  length;  oi  =  circular  wave  frequency 
z'  -  vertical  coordinate,  from  undisturbed  water 
surface,  positive  downwards 
x'  =  axis  fixed  in  space 
t  =  time. 

The  x'  coordinate  of  a  point  in  the  x-y  plane  can  be  defined  by: 

x'  =  -(x+Vt)  cos  8  +  y  sin  8  (2) 

The  surface  wave  elevation  n  (positive  upwards)  can  be 
expressed  as  follows: 


direction  of  ship  travel 
at  speed,  V 


Fig.  1.  Wave  and  Ship  Axes  Convention 
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n  =  — 


i.  ~  w 

9  l  at  i 


=  a  sin  k  (x'  +  ct) 


(3) 


z  '=0 


since 


<=2  =  ; 

where  g  =  acceleration  of  gravity. 

In  x-y  coordinates,  relative  to  the  ship,  we  have: 


n  =  a  sin  k  [-x  cosB  +  y  sinS  +  (c-V  cosB)t],  (4) 


Dn 

Dt 


i_  -v  — 
TE  3x 


n  (x,t) 


and  n  =  §£ 


=  akc  cos  k  [-x  cosB  +  y  sinB  +  (c-V  cosB)t] (5) 
=  -akg  sin  k  [-x  cosB  +  y  sinB  +  (c-V  cosB)t](6) 


The  results  of  the  equations  of  motion  and  the  wave- induced 
moments  will  be  referenced  to  the  wave  elevation  n  at  the  origin 
of  the  x-y  axes,  which  is 


n  **  a  sin  k  (c-V  cosB)  t 


(7) 


or 


n  =  a  sin  w  t 
e 


where 


2  TT 


we  =  (c-V  cosB) 


(8) 


and  ue  is  known  as  the  circular  frequency  of  encounter.  The  quantity 
ue  is  generally  positive,  and  only  for  following  waves  (90°<B<90°  ), 


where  the  ship  is  overtaking  the  waves,  is  a>e  negative. 


Vertical  Plane  Equations 


The  coupled  equations  of  motion  in  the  vertical  plane  for  heave, 
z  (positive  downwards),  and  pitch,  6  (positive  bow-up) ,  in  keeping 
with  the  revised  axes  convention,  are  given  as: 


(9) 


h 

A 


6 


••  A7 

V  -  -  a§  x  +  Mw 


where  m  =  mass  of  ship 

I  =  mass  moment  of  inertia  of  ship  about  y  axis 

^  =  local  sectional  vertical  hydrodynamic  and  hydrostatic 
force  on  ship 

x  ,x.  =  coordinates  of  stern  and  bow  ends  of  ship, 
s  respectively 

Zw,Mw  =  wave  excitation  force  and  moment  on  ship. 

The  general  hydrodynamic  and  hydrostatic  force  is  taken  to  be: 

g|  =  -  (Z-X8+V6 )J -N^ (z-x* +V0) -pgB* (z-x6 )  (] 

where  p  =  density  of  water 

A'33  =  local  sectional  vertical  added  mass 

N'z  =  local  sectional  vertical  damping  force  coefficient 
B*  =  local  waterline  beam 
_2  -3 

and  N'  =  pg2A  |u  |  (] 

2  6 

with  A  =  ratio  of  generated  wave  to  heave  amplitude 

for  vertical  motion- induced  wave. 

Values  of  the  two-dimensional  A33  and  A  terms  are  calculated  by 
the  method  of  Grim  [7]  for  the  equivalent  Lewis  forms  at  each 
section. 

Expanding  the  derivative  in  Eq.  (11),  we  obtain: 

=  -  A'3  (Z-X0+2V0)  -  (S-X0+V0) 

-  pgB*(z-x0)  (- 

The  equations  of  motion,  (9)  and  (10)  are  transformed  into  the 
familiar  form  as  follows: 


a'z  +  bz  +  c*z  -  de  -  ee  -  g'0  =  Z 


(14) 
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where  h  =  mean  section  draft.  Substituting  the  expressions  for  n, 
n  and  n  from  Eqs.  (4),  (5)  and  (6),  with  y=0,  and  incorporating  an 
approximate  factor  for  short  wave  lengths  (by  carrying  out  an 
integration  over  the  lateral  extent  of  the  ship,  in  terms  of  the 

n*  n* 

y-coordinate  extending  from  -j-  to  ) ,  leads  to 


where  the  latter  factor  in  Eg.  (20)  represents  this  short  wave 
length  factor.  The  value  of  h  is  approximated  by: 

h  =  HCg  (21) 

where  H  =  local  section  draft 

C  =  local  section  area  coefficient, 
s 

The  steady-state  solution  of  the  equations  of  motion,  at  each 
particular  regular  wave  length,  is  obtained  by  conventional  methods 
for  second  order  ordinary  differential  equations  (using  complex 
notation).  The  solutions  are  expressed  as: 


z  =  z  sin(w  t+6 ) 
o  e 

(22) 

0  -  dQ  sin(uet+e) 


where  the  zero  subscripted  quantities  are  the  motion  response 
amplitudes  and  <$,  e  are  the  phase  angle  differences,  i.e.  leads 
with  respect  to  the  wave  elevation  in  Eq.  (7) , 

Having  obtained  solutions  for  the  motions  in  the  vertical 
plane,  the  wave-induced  vertical  bending  moment  can  then  be 
calculated.  The  bending  moment  is  found  from  the  total  loading 
at  each  section.  This  is  made  up  of  the  loads  due  to  inertia 
(ship  mass),  hydrodynamic  and  hydrostatic  forces,  and  the  direct 
wave  loads.  The  total  local  vertical  loading  is  then  given  by: 
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dfz  ••  ••  dZ  dZw 

33T  =  “  6m  (2_xe>  +  a x  +  33T 


where  6m  =  local  mass. 

Eq.  (23)  is  simply  the  summation  of  inertial,  hydrodynamic,  hydro¬ 
static  and  wave  excitation  forces.  The  latter  terms  are  given  in 
Eqs.  (13)  and  (20).  The  vertical  bending  moment  at  any  location 
xq  along  the  ship  length  is  then  given  by: 


BW  - 


°r  '“-V  ffiT  a* 


and  is  expressed  in  a  form  similar  to  the  motions,  i.e. 


BM  =  BM„  sin  (ui  t+o) 
z  zo  e 


Lateral  Plane  Equations 

The  coupled  equations  of  motion  in  the  lateral  plane  for  sway, 
y  (positive  to  starboard),  yaw,  i|/  (positive  bow- starboard)  ,  and  roll, 
6  (positive  starboard-down),  are  given  as: 


my  = 


a?  dx+Yw 


I  li  -I  6  = 

zv  xzv 


a!  xdx+Nw 


-I 

xz 


dx-mg  GM  <f>+Kw 


where  lz  =  mass  moment  of  inertia  of  ship  about  z  axis 
Ix  =  mass  moment  of  inertia  of  ship  about  x  axis 
I  =  mass  product  of  inertia  of  ship  in  x-z  plane 

XZ 


J  In,  hJ-UK.tU.  • 
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=  local  sectional  lateral  hydrodynamic  force  on  ship 

jpy  =  local  sectional  hydrodynamic  rolling  moment  on  ship 

Y  ,  N  ,  K  =  wave  excitation  force  and  moments  on  ship 
w  w  _w 

GM  =  initial  metacentric  height  o*  ship  (hydrostatic) . 

The  cross  inertial  terms,  involving  I  ,  the  mass  product  of  inertia, 
are  usually  small  but  necessary  for  tfte  equilibrium  balance  of  forces 
and  moments.  The  hydrodynamic  force  and  moment  in  the  above 
equations  are  given  by: 


[Ms  (y+xJ>-V«|<)  -Frs$J  -Ng  (y+xji-' 


Vip)  +  Nrs4» 


+  <5G  (Mg$)  +  OG  Ng$ 


(y+xifi-Vii»)J  -Nr$+Ng^  (y+x*-V<|>) 


-  0G  §E  ‘  ^  Ns^  "  ™ 


where  OG  =  distance  of  ship  C.G.  from  waterline,  positive  up 

M  =  sectional  lateral  added  mass 
s 

N  =  sectional  lateral  damping  force  coefficient 
s 

M  =  sectional  added  mass  moment  of  inertia  due  to  lateral 
motion 

N  =  sectional  damping  moment  coefficient  due  to  lateral 
motion 

I  =  sectional  added  mass  moment  of  inertia 
Nr  =  sectional  damping  moment  coefficient 

F  =  sectional  lateral  added  mass  due  to  roll  motion 
rs 

N  =  sectional  lateral  damping  force  coefficient  due  to 
rs  roll  motion 

and  the  sectional  added  mass  moments  and  damping  moment  coeffi¬ 
cients  are  taken  with  respect  to  an  axis  at  the  waterline. 

Values  of  these  sectional  hydrodynamic  properties  for  the 
equivalent  Lewis  form  at  each  section,  as  functions  of  the 
frequency  of  oscillation,  can  be  calculated  by  the  method  of 
Tasai  (8)  based  on  the  potential  theory  solution.  It  has  been 
shown  by  Vugts  [9]  that  such  potential  theory  results  for  the 
lateral  and  rolling  modes,  which  ignore  viscous  and  surface 
tension  effects,  are  in  good  agreement  with  experimental  results 
except  for  the  roll  damping  moment.  In  addition,  the  influence  of 
bilge  keels,  which  are  usually  used  but  not  considered  up  to  this 
point,  is  expected  to  be  primarily  upon  the  roll  damping  moment. 
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In  order  to  account  for  the  above  effects,  that  is  the  viscous 
effect  and  the  bilge  keel  effect  upon  the  roll  damping  moment,  an 
adjustment  is  made  to  the  potential  theory  result.  Roll  motion  is 
generally,  for  conventional  merchant  hull  forms,  a  very  lightly 
damped  response.  This  means  that  at  resonance,  i.e.  at  the  natural 
roll  response  frequency,  the  damping  value  is  important  in  limiting 
large  roll  responses,  but  that  at  frequencies  away  from  resonance 
the  amount  of  damping  hardly  affects  the  roll  response  at  all.  Thus 
it  is  most  important  to  determine  the  proper  value  of  the  roll 
damping  moment  at  the  resonant  frequency,  while  at  other  frequencies 
away  from  resonance  its  influence  is  almost  negligible.  The 
adjustment,  or  addition,  to  the  roll  damping  moment  is  made  so  that 
at  the  resonant  frequency  the  total  roll  damping  is  a  particular 
fraction  of  the  critical  roll  damping.  This  fraction  is  estimated, 
or  known  by  experimentation,  to  produce  the  proper  roll  response 
at  resonance.  This  approach  was  employed  by  Vugts  [10]  and  verified 
experimentally  for  the  rolling  motions  of  a  cylinder  of  rectangular 
cross-section  in  regular  beam  waves.  Therefore,  we  have: 


NJ  -  SCC/L  -  V"*>  <311 

where  N*  =  sectional  damping  moment  coefficient  due  to 

r  viscous  and  bilge  keel  effects 

=  fraction  of  critical  roll  damping  (empirical  data) 

Cc  =  critical  roll  damping 

L  =  ship  length  (L  =  xb~xg) 

=  natural  roll  (resonant)  frequency 

N  (w.)  =  value  of  N  at  frequency  of  w.. 
r  <p  r  9 

This  procedure  is  still  linear,  with  the  empirical  value  of  the 
damping  at  resonance  reflecting  an  average  or  equivalent  linear 
value  that  can  be  applied  in  an  approximate  manner.  Since  the  main 
concern  of  this  study  is  determining  structural  loads,  and  the 
influence  of  roll  motion  per  se  must  be  explored  in  the  investi¬ 
gation  itself,  the  use  of  this  method  of  representation  is 
considered  sufficiently  valid  for  this  purpose.  The  critical 
roll  damping  can  be  expressed  in  terms  of  the  natural  roll 
frequency  as  follows: 


with 


(32) 


where  the  integral  is  over  the  ship  length 
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Expanding  the  derivatives  in  Eqs.  (29)  and  (30),  and  including 
the  above  additional  roll  damping  moment,  we  obtain: 

dY  ...../  dM _  \  .  _  \  .. 

=  "  Ms(y+x^-2V^)  +  V  -Nj  (y+x^-V*)  +  Frs+  OG  Mj  4> 

_  /dF  _  dM  \  1  . 

+  |Nrs  +  OG  Ns-  V  +  OG  g/j  I  ♦  (33) 

dK  r  _  /  \  ~ I  r  /  dl  dM 

3x  =  ”  JCr$+  0G  |Ms*+  Frs+  0G  Msj  ♦  +  v  SsT  +  0G  "cRT/ 


dl  _  dM 

g—  +  °C  _g_i 


-  0G  Nrs+  Ns*+  0G  Ns  +  0G  v  "ST  +  0G  35r, 


-Nr-  N*  4 


+  Mg(ji+  OG  Ms  (y+Xt}l-2V^)  +  N_  +  OG  Ns-V  +  UG 


(y+xJr-ViJ)) 

The  equations  of  motion,  (26) ,  (27)  and  (28)  are  transformed 
into  this  familiar  form: 


all^+a12^+a14^+a15^+a16^+a17^+a18^  =  Vw 
a21^+a22^+a24^+a25^+a26^+a27^+a28^  =  Nw 
a31^+a32^+a34^+a35^+a36^+a37^+a38^+a39^  =  Kw 


The  coefficients  on  the  left-hand  sides  are  then  defined  by: 


aH  =  m  +  |Msdx  ,  a12  =  J^dx-V  |d(Ms)  , 

a14  =  jMsxdx  '  ais  =  |Ngxdx-2V  |Msdx-V  |xd(Mg) 

a16  =  -Va12  '  a17  =  -  fFrsdx  "  6G  |Msdx  ' 

a18  =  "  fNrsdx  +  |d(Ms)-OG  |Nsdx  +  V  jd(Frg) 


a21  -  fMsxax 
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,  a^2  =  |nsxcIx-V  jxd  (Ms)  , 

=  Iz+  fMsx2dx  ,  a25  =  jNgx2dx-2V  |Msxdx-V  Jx2d(Ms), 

_ Vs  a  =  -I  -If  xdx  -OG  M  xdx  , 

-  Va22  '  27  xz  Jr s  J  s 

=  -  |Nrsxdx+OGv|xd(Ms)  -  OG  j^gXdx+V  jxd(Frg)  . 

-  -  -  53  JMsd*  ' 

.  -  Jti^dx  -  OG  Jtjgdx  +V  Jd(Mg^)  *  V  OG  Jd(M,.)  , 


1  Adx  -  OG  M  dx  , 

IS 


a34  -  -Zxz  -  K»XdX  "  °G  iM=XdX  ' 


=  -  jllg^xdx  -  OG  |Ngxdx  +V  |xd(Mg  +  )+V  OG  jxd  (Mg) -2Va31  , 


a36  =  ”Va32' 


=  Ix+  jlrdx  +  OG  jMg^dx  +  OG  |Frsdx  +  OG2  jt^dx  , 
=  |  (Nr+N*)dx  +  OG  |Ns<j)dx  +  OG  Jl^dx  +  OG  2  JNgdx 

-  V  Jd(Ir)+OG  Jd{MS(f)+OG  jd(Frs)+OG2  Jd(Mg) 


a39  =  mg  GM 


where  all  the  indicated  integrations  are  over  the  ship  length. 

The  wave  excitation ,  the  right-hand  sides  of  Eqs.  (35) t  is 
given  by: 


Yw  = 


(39) 
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(44) 
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The  steady-state  solution  of  the  lateral  plane  equations  of 
motion,  at  each  particular  regular  wave  length,  are  expressed  as: 


y  =  yQ  sin  (wet  +  *)  (48) 


iji  =  i|>o  sin  (w^t  +  a)  (49) 


=  4>o  sin  (uet  +  v)  (50) 


where  the  zero-subscripted  quantities  are  the  motion  response 
amplitudes  and  k,  a  and  v  are  phase  angle  leads  with  respect  to 
the  wave  elevation. 

The  local  lateral  (force)  and  rotational  (moment)  loadings 
derived  in  a  manner  similar  to  the  vertical  loading,  are  given  by: 


(51) 


where  ?  =  local  center  of  gravity  (relative  to  ship  C.G.), 
positive  down 

y  =  local  mass  gyradius  in  roll 

and  the  hydrodynamic  and  wave  excitation  terms  are  given  in  Eqs. 
(33),  (34),  (46),  and  (47).  While  the  local  lateral  loading  is 
directly  comparable  to  the  local  vertical  loading,  including 
inertial,  hydrodynamic  and  wave  excitation  forces,  the  local 
rotational,  or  torsional,  loading  must  in  addition  account  for  the 
static  rotational  momen£,  due  frh,e  ^aiiiteJ^^etaeeftttfic^heicffil?-* 

♦  -teft-SlT  difa* local  rse*ctional)  basis. 

Finally,  the  wave-induced  lateral  bending  moment  and 
torsional  moment  at  any  location  x  along  the  ship  length  are  then 
given  by: 


I'Swm.-vrr--'.*  ■'■'■RMn.'vw  >w,iyiiav«flnv<u« 


I 


BMy(X0)  = 


df 

(54"xo)  a/  dx 


TMx(xo)  = 


and  again  they  are  expressed  in  this  form: 


BMy  =  BMyo  sin  (“et  +  T) 


™x  =  ™xo  sin  (“et  +  v) 


The  parameters  defining  the  ship  mass  distribution  must  meet 
certain  constraints.  The  requirement  on  ?,  the  local  vertical  mass 
center,  is: 


6m?  dx  =  0 


since  ?  is  measured  relative  to  the  ship  C.G.,  and  all  first  moments 
about  that  point  must  sum  to  zero,  by  definition.  Similarly,  the 
requirement  on  v,  the  local  roll  gyradius,  is: 


r  ? 


The  product  of  inertia  in  the  x-z  plane  is  then  defined  by: 


6mx?dx 
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We  should  note  here  that  it  is  usual  practice  in  model  test  work 
that  each  overall  segment,  or  portion,  of  the  model  is  ballasted 
to  the  same  overall  specified  V.C.G.  and  roll  ryradius.  However, 
data  concerning  the  variation  of  r,  and  y  with  length  is  usually 
not  available. 


Irregular  Sea  Equations 


All  of  the  results  obtained  in  the  preceedinq  analyses  have 
been  appropriate  to  conditions  of  regular  sinusoidal  unidirectional 
waves,  which  occur  only  in  model  test  tanks.  In  a  realistic  seaway, 
waves  appear  randomly,  and  the  motions  and  structural  responses 
of  a  ship  in  such  waves  also  have  a  random  nature.  In  order  to 
characterize  the  random  ship  responses,  the  energy  spectra  of  the 
responses  are  employed.  Each  spectrum  is  a  measure  of  the 
variation  of  the  squares  of  the  amplitudes  of  the  various  sinusoidal 
components  of  the  particular  random  response,  presented  as  a  function 
of  .the  wave  frequency.  The  spectral  technique  for  analyzing  random 
irregular  time  histories  of  motion  and  structural  response  is 
applicable  to  linear  systems  only,  since  in  that  case  a  unique 
response  amplitude  operator  is  obtained.  The  spectral  techniques 
evolve  as  a  result  of  linear  superposition,  as  originally  developed 
in  [11] ,  of  the  responses  to  individual  frequency  components 
contained  in  the  wave  excitation. 


The  surface  wave  system,  which  is  defined  by  the  wave  energy 
spectrum,  is  considered  to  be  a  separable  function  of  wave  frequency 
and  direction,  with  limits,  as  follows: 


S  (u,w)  =  S 2  ( p )  for  0<u><»  (59) 

.  IT  71 

and  ~  2  —  v  -  2 

where  S  (w ,u)  =  directional  spectrum  of  the  seaway  (short 
crested  sea  spectrum) 
to  =  circular  wave  frequency 

u  =  wave  direction  .relative  to  predominent  direction 
S^(w)  =  frequency  spectrum  (long  crested  sea  spectrum) 

S2(u)  =  spreading  function. 

The  mean  squared  wave  amplitude  is  a  basic  measure  of  the  total 
energy,  or  intensity,  of  the  particular  sea  spectrum.  It  is 
obtained  simply  as  the  integral  of  all  the  various  components,  in 
continuous  form,  as: 

'  »  IT 

—  2 

,  a2  =  S(w,u)  du  du  (60) 

TT 


0 


'  “  *LL  V1  44  J 
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where  a2  =  mean  squared  wave  amplitude,  or  variance  of  the  wave 
time-history  record.  Since  the  spreading  function  depends  on 
relative  wave  direction  only,  it  is  usual  to  impose  the  following 
constraint: 

rir 

2 

S2(u)  du  =  1.0  (61) 

71 

;  Therefore,  we  can  define  the  mean  squared  wave  amplitude  in  terms 

i  of  the  long  crested  sea  spectrum  as: 


.2  - 


(w)  do) 


(62) 


Other  statistical  parameters  of  interest  for  the  sea  spectrum,  and 
similarly  for  any  response  spectrum,  can  be  obtained  from  the  mean 
squared  amplitude,  or  variance,  of  the  particular  random  variable. 
For  the  waves,  we  have: 


arms 
aavg 
al/3  = 
al/10= 


1/2 

(a2) 

(63) 

1.25  a__ 

(64) 

rms 

2.0  a„mc. 

(65) 

rms 

2.55  a 

rms 

(66) 

where  arms  =  root“mean”st3uarec*  wave  amplitude 

a_„_  =  average  (statistical)  wave  amplitude 

ctvy 

a./3  =  significant  (average  of  1/3  highest) 

'  wave  amplitude 

a-jy^0=-  average  of  1/10  highest  wave  amplitude. 

Various  long  crested,  or  unidirectional,  sea  spectra  have  been 
proposed  over  the  years  as  representative  of  realistic  conditions 
at  sea.  Three  spectral  formulations  in  popular  usage  among  various 
investigators  in  the  field  are  given  below,  for  reference. 
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Neumann  Spectrum  (1953) :  This  frequency  spectrum  [12]  can 
be  specified  by: 

S^w)  =  0.000827  g2ir3u“6e“2g2u“2u"2  (67) 

where  U  =  wind  speed.  The  constant  given  here  is  one  half  that 
originally  specified  by  Neumann,  so  that  this  spectrum  satisfies 
Eq.  (62).  Thus,  originally  the  Neumann  spectrum  required  only  a 
factor  of  /2~  in  Eq.  (65),  instead  of  2.0. 

Pierson-Moskowitz  Spectrum  (1964)  ;  This  is  given  [13]  by: 

S1(lo)  =  0.0081  g2w-5e''*74g4“)~4u”4  (68) 

and  was  derived  on  the  basis  of  fully  arisen  seas. 

Two  Parameter  Spectrum  (1967):  This  spectrum  is  intended  for 
use  in  conjunction  with  ^observed'1  wave  height  and  period,  which 
are  then  taken  to  be  the  significant  height  and  mean  period.  This 
spectrum  is  similar  to  that  adopted  by  the  I.S.S.C.  (1967)  [14]  as 
"nominal",  except  that  it  is  expressed  here  in  circular  wave 
frequency  instead  of  frequency  in  cycles  per  second: 


S^(uj)  =  A’Bu“5e  — w 


where  A  =  0.25  (*^3)  2 

B  »  (0.817  l^-)4 
T 

H1/3  =  significant  wave  height  (=2. °a1^3) 

A» 

T  =  mean  wave  period. 

The  spreading  function  can  be  expressed  for  long  crested, 
or  unidirectional,  seas  as  follows: 

S2(u)  =  6(p)  (70) 

where  5 (  )  =  delta  function. 

For  short  crested  seas,  various  spreading  functions  have  been 
suggested  and  developed.  Perhaps  the  most  widely  used,  and  a 
compromise  among  the  proposed  forms,  is  the  cosine-squared 
spreading,  expressed  as: 

2 

S- (u)  =  —  cos2m 


(71) 
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Having  defined  the  wave  excitation,  or  sea  spectrum,  in  the 
forms  as  given  above,  the  energy  spectrum  of  the  motions  or 
structural  responses  can  be  calculated.  In  line  with  the  linear 
assumption  for  all  responses,  and  employing  the  principle  of  wave 
superposition,  a  response  spectrum  is  obtained  by: 


(w,m) 


( U) ,  U  ) 


(to,p) 


(72) 


where  S^(u>,u)  =  response  spectrum,  for  a  particular  response 

T. (w,u)  =  response  amplitude  operator  (amplitude  of  i-th 
1  response  per  unit  wave  amplitude) . 

We  then  have,  similar  to  the  wave  amplitude: 


oo 

(u) 

0 


where  a^2  =  mean  squared  response  amplitude.  The  term  in  square 
brackets  in  the  integral  above  is  the  mean  squared  response 
amplitude  for  long  crested  seas  at  a  particular  heading  u, 
relative  to  the  predominant  wave  direction.  The  other  statistical 
parameters  of  interest  for  each  response  can  be  obtained  from  the 
mean  squared  amplitude  by  use  of  equations  similar  to  Eqs.  (63) 
to  (66). 


EXPERIMENTAL  DATA 

In  order  to  evaluate  the  analytical  methods  presented  for 
the  calculation  of  wave-induced  moments,  the  ref  "ilts  of  such  cal¬ 
culations  are  to  be  compared  with  experimental  results.  Experi¬ 
mentation  on  ship  models,  under  controlled  laboratory  conditions, 
for  the  determination  of  wave-induced  moments  is  a  relatively  new 
procedure.  Lewis  [15]  first  presented  such  results  for  vertical 
bending  in  1954.  These  initial  tests  were  limited  to  head  and 
following  seas  directions.  The  tests  measured  midship  bending 
moments  only. 
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Subsequently,  the  experimental  procedures  were  expanded  to 
cover  a  greater  degree  and  range  of  relevent  parameters,  approach¬ 
ing  description  of  responses  under  various  realistic  conditions. 
Among  the  earliest  tests  conducted  at  oblique  wave  headings  in 
order  to  yield  both  vertical  and  lateral  bending  moments,  and 
torsional  moments  as  well,  were  those  of  Numata  [16]  conducted 
at  Davidson  Laboratory  on  a  T-2  tanker  model.  This  work  was 
quickly  followed  by  an  extensive  series  of  tests  on  Series  60 
models  by  Vossers,  et  al,  [17]  at  the  Netherlands  Ship  Model  Basin 
(NSMB) ,  reported  in  1961.  In  addition,  tests  also  have  been  con¬ 
ducted  for  determination  of  wave-induced  loads  at  points  along  the 
hull  other  than  midships.  However,  at  this  time  the  total  amount 
of  all  such  data  is  not  very  large  and  some  experimental  problems 
still  exist.  Very  little  data  has  been  collected  with  regard  to 
torsional  moments,  and  therefore  the  emphasis  in  the  comparison 
to  follow  will  be  upon  the  vertical  and  lateral  bending  moments 
at  midships. 

A  fairly  intensive  test  series  was  reported  by  Wahab  [18]  in 
1967.  These  tests  of  only  one  Series  60  hull  form,  with  block 
coefficient  of  0.80,  were  conducted  over  a  large  and  exhaustive 
range  of  regular  wave  lengths  and  wave  angles.  Measurements  in¬ 
cluded  vertical  and  lateral  bending  moments,  plus  vertical  and 
lateral  shears  and  torsional  moment,  all  at  midships.  Recently 
the  Ship  Structure  Committee  has  supported  additional  experimental 
work  at  Davidson  Laboratory  that  is  related  to  other  full-scale 
measurement  projects.  The  model  tests  have  been  reported  by 
Chiocco  and  Numata  [19]  for  the  "Wolverine  State,"  and  by  Numata 
and  Yonkers  [20]  for  the  Mariner-class  "California  Bear." 

With  regard  to  the  comparison  between  such  experimental  data 
and  the  projected  calculations,  certain  conditions  of  the  model 
tests  should  be  recognized.  The  bulk  of  the  test  data  to  be  used 
in  this  comparison  are  the  results  of  model  tests  in  regular  waves 
at  oblique  headings,  referenced  above.  Such  tests  are  conducted 
by  using  a  fairly  free-running  self-propelled  model.  The  model 
must  then  have  an  operational  rudder  which  is  used  to  maintain 
the  model  along  the  prescribed  wave-to-course  angle.  In  more  re¬ 
cent  tests  of  this  type,  the  rudder  is  controlled  by  an  automatic 
procedure  and/or  device  based  on  yaw  and  sway  motions  that  are 
sensed  by  elements  on  the  model,  while  in  some  early  tests  the 
rudder  was  controlled  manually.  In  regular  wave  tests,  it  would 
appear  that  the  rudder  movements  could  contain  significant  encounter 
frequency  content.  In  any  event  it  is  clear  that  the  rudder 
action  influences  the  model  motion  responses  under  such  conditions. 
Furthermore,  the  rudder  forces  generated  in  this  manner  contribute 
directly  to  the  total  loading  distribution  on  the  hull,  which  is 
assumed  to  be  in  equilibrium.  Since  the  lever  arm  of  the  rudder 
forces  is  large  for  moments  at  midships,  it  appears  that  rudder 
forces  can  significantly  affect  the  lateral  bending  and  torsional 
moments.  To  the  extent  that  the  use  of  the  rudder  affects  the 
overall  ship  motion  response  in  oblique  seas,  the  vertical  bending 
moment  also  can  be  influenced,  but  to  a  much  smaller  degree.  The 
calculations,  based  on  the  analytical  method  presented  earlier, 
do  not  include  any  rudder  force  and  moment  effects. 
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The  above  discussion  of  rudder  effects  only  points  out  a 
difference  between  experimental  conditions  and  the  proposed 
calculations.  Another  point,  and  one  of  perhaps  equal  importance 
but  not  directly  bearing  on  the  subject  comparison  at  hand,  is 
whether  such  model  scale  rudder  forces  and  control  techniques  are 
representative  of  full-scale  effects.  Questions  of  scale  effect 
and  response  times  enter  into  this  problem,  and  will  not  be 
considered  here.  The  point  is,  however,  that  the  calculations 
are  being  compared  with  experiments  which  include  additional  un¬ 
accounted  effects,  which  are  not  necessarily  realistic  with 
regard  to  full-scale  behavior. 

Another  aspect  of  the  experimental  conditions  also  is 
significant  with  regard  to  the  comparison  with  calculations.  In 
the  experiments  at  oblique  wave  angles ,  it  is  noted  that  the 
model's  mean  heading  angle  differs  from  the  mean  wave-to-course 
angle,  the  difference  being  referred  to  as  the  leeway  angle.  The 
leeway  angle  appears  to  be  due  to  the  non-zero  mean  lateral  forces 
and  moments  imposed  by  the  waves.  It  is  greatest  at  low  speeds 
in  relatively  short  wavelengths.  Thus,  for  example,  in  an  exper¬ 
iment  at  a  wave-to-course  angle  of  120°,  bow  seas,  the  actual 
average  heading  of  the  hull  to  the  waves  may  be  as  high  as  135° 
[16].  The  analytical  methods  take  no  account  of  such  mean,  or 
drift,  forces  and  moments,  so  that  in  the  calculations  the  leeway 
angle  is  assumed  to  be  zero.  Since  no  account  of  the  leeway  angle 
is  made  in  the  computations,  and  the  wave-to-course  angle  used  in 
the  computations  is  the  nominal  value  prescribed  in  the  tests,  the 
influence  of  the  actual  heading  of  the  ship  relative  to  the  waves 
is  not  accounted  for  properly.  In  the  model  tests,  the  mean 
wave  forces  and  moments  which  cause  the  leeway  angle,  and  the  mean 
hydrodynamic  forces  and  moments  resulting  thereby,  are  supposedly 
in  balance  with  the  force  and  moment  from  a  non-zero  mean  rudder 
angle.  That  is,  it  is  usually  necessary  to  apply  a  mean  rudder 
angle  in  order  to  keep  the  model  on  a  prescribed  mean  course,  but 
with  a  particular  resulting  leeway  angle. 

Obviously,  these  forces  and  moments  have  some  effect  upon 
the  motion  responses  of  the  model  and  therefore  upon  the  measured 
moments.  The  extent  and  nature  of  such  effects  are  unknown,  al¬ 
though  the  only  important  effects  will  be  those  forces  at  the 
frequency  of  encounter  in  the  regular  wave  tests.  However,  in 
the  reports  of  the  experimental  work  little  or  no  significance 
is  given  to  these  forces.  The  details  of  the  rudder  and  control 
system  are  not  described.  Rudder  motion  is  not  given,  and  even 
leeway  angle  is  not  always  reported.  Thus,  at  this  point  in  the 
development,  the  experimental  inputs  for  comparison  with  a  full 
analytical  treatment  of  rudder  forces  and  mean  wave  forces  and 
moments,  if  such  were  desired,  are  not  yet  available.  The  effects 
then  of  leeway  angle  and  rudder  forces  may  turn  out  to  be  small 
in  many  cases,  but  they  must  still  be  recognized  as  an  unknown 
element  in  the  comparison. 
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COMPARISON  OF  CALCULATIONS  WITH  EXPERIMENT 

All  of  the  calculations  of  wave-induced  moments  were  done 
by  use  of  a  dig’  .1  computer  program  developed  in  the  course  of 
this  work,  and  fully  described  elsewhere  [3] .  The  program  follows 
the  analytical  methods  presented  in  this  report  and  its  pre¬ 
decessors  [1,  2] .  The  calculations  of  the  midship  wave-induced 
moments  were  carried  out  for  hull  forms,  mass  distributions  and 
test  conditions  corresponding  to  the  bulk  of  the  experimental 
data  cited  previously.  In  general,  sufficient  data  was  available 
in  the  model  test  reports  with  regard  to  the  full  description  of 
the  necessary  significant  parameters  for  input  to  the  computer 
calculation.  However,  as  pointed  out  previously,  no  data  was 
available  with  regard  to  the  longitudinal  distribution  of  c,  the 
local  vertical  center  of  gravity,  and  y,  the  local  roll  gyradius. 
These  parameters  can  be  expected  to  affect  the  lateral  bending 
moment  in  the  region  of  roll  resonance  only,  and  also  for  the 
torsional  moment.  In  some  cases  a  reasonable  approximation  to 
the  vertical  center  of  gravity  distribution  was  employed,  corres¬ 
ponding  to  the  usual  model  test  ballasting  methods.  In  these 
cases,  the  lateral  bending  moment  calculation  results  were  seen 
(via  numerical  tests),  to  be  sensitive  to  this  distribution  in  the 
region  of  roll  resonance.  The  use  of  a  reasonable  approximation 
generally  yielded  results  which  were  in  better  agreement  with  the 
experimental  results. 

In  order  to  simplify  the  presentation  of  the  results  of  the 
computations,  and  comparison  with  model  test  data.  Table  1  has 
been  prepared.  It  lists  the  calculations  to  be  presented  herein, 
together  with  the  reference  for  the  experimental  results.  For 
each  of  the  five  sets  of  calculations.  Tables  2  to  6  give  the 
basic  hull  form  and  mass  distribution  data  used,  based  on  the  input 
values  specified  and  inherent  assumptions  in  the  computer  program. 
Also  shown  in  Table  1  are  the  roll  damping  fractions  used  in  the 
computations  for  lateral  plane  motions ,  and  the  figure  numbers 
which  give  the  results,  including  comparison  with  corresponding 
experimental  data. 

Primarily  the  comparison  is  made  for  the  Wolverine  State  data 
[19]  and  the  Series  60,  block  0.80  hull  data  of  Wahab  [18].  These 
represent  more  recent  tests  of  this  type,  where  experimental  pro¬ 
cedures  are  perhaps  more  refined  compared  to  earlier  tests.  The 
Wolverine  State  comparison  is  for  two  different  hull  loadings,  two 
speeds,  and  over  a  fairly  wide  range  of  wave  angles  and  wavelengths. 
The  Series  60,  block  0.80  hull  comparison  is  at  one  loading  and 
speed,  but  the  experimental  data  cover  a  wide  range  of  wave  lengths 
and  angles  more  intensively.  The  comparison  is  also  presented  for 
the  Series  60,  block  0.70  hull  data  (NSMB,  1961)  and  the  T-2  Tanker 
Model  (Davidson  Lab.,  1960)  so  that  a  wider  range  of  hull  forms 
and  test  conditions  can  be  covered.  From  Tables  2  to  6,  it  can  be 
seen  that  twenty  stations  along  the  ship's  length  were  generally 
used  to  define  the  hull  form  and  mass  distribution.  This  is  con¬ 
sidered  an  appropriate  number,  compared  with  other  numerical  aspects 
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Table  1.  Calculations  Reference  Data 


Model  Description 

Model  Test 
Reference 

Basic 

Data 

S 

Calculation 

Results 

1 

Wolverine  State, 
Full  Load 

[191 

Table 

2 

0.05 

Figs.  2,  3 

2 

Wolverine  State, 
Light  Load 

[19] 

Table 

3 

0.05 

Figs.  2,  3 

3 

Series  60, 

Block  .80  Hull 

[18] 

Table 

4 

0.10 

Figs.  4-7 

4 

Series  60, 

Block  .70  Hull 

[171 

Table 

5 

0.05 

Figs.  8-11 

li 

T-2  Tanker  Model 

[161 

Table 

6 

0.05 

Figs.  12 

Table  2.  Basic  Data  for  WOLVERINE  STATE,  Full  Load 
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GRAVITY  «  32, 

,175000 

STATION  ream 

ACE*  COEE 

.  OOAET 

Z-bar 

HEIGHT 

ZETA 

gyr.roll 

n.oo 

0.0000 

0.0000 

o.oooo 

0.0000 

104,0000 

-21,0000 

26.6000 

1.00 

11.2000 

.8060 

30.0000 

13.0600 

306,0000 

-20,6346 

26.6000 

2.00 

24,4000 

.84R0 

30.3000 

13.4800 

532.0000 

-20.0000 

26,6000 

3.00 

38.6500 

.8840 

30,0000 

13,8400 

862,0000 

-17.0000 

26.6000 

*.00 

SI. 4000 

.0130 

30.0000 

14,1300 

1060.0000 

-12.0000 

26.6000 

5,00 

61,3000 

.0330 

.30,0000 

14,3300 

1201.0000 

-6,0000 

26.6000 

0.00 

67.4000 

.0570 

30.0000 

14.5700 

1310.0000 

2,0000 

26.6000 

T.00 

70.5000 

.0750 

30.0000 

14,7500 

1399,0000 

8.0000 

26.6000 

B.OO 

71,5000 

.0840 

30.0000 

14,8400 

1418.0000 

10.0000 

26.6000 

9.00 

71.5000 

.0890 

30.0000 

14.8900 

1428.0000 

12.0000 

26.6000 

10.00 

71 .5000 

.0890 

30,0000 

14,8900 

1442.0000 

12.0000 

26.6000 

11.00 

71.5000 

.0840 

30.0000 

14,8400 

1446.0000 

12.0000 

26.6000 

12.00 

71.5000 

.0670 

30.0000 

14.6700 

1395.0000 

10.0000 

26.6000 

13.00 

71.5000 

.0310 

30.0000 

14,3100 

1296,0000 

8.0000 

26.6000 

1*.00 

71.2000 

.8630 

30.0000 

13.6300 

1079.0000 

2.0000 

26,6000 

15.00 

69,4000 

.7720 

30,0000 

12,7200 

791,0000 

-6.0000 

26,6000 

16.00 

63.4000 

.6730 

30,0000 

11.7300 

716.0000 

-12.0000 

26.6000 

17.00 

52.3000 

.5830 

30.0000 

10.8300 

772.0000 

-17.0000 

26,6000 

18.00 

37.4000 

.4950 

30.0000 

9.9500 

593.0000 

-19.0000 

26,6000 

19.00 

21.9500 

.3860 

30.0000 

8,8600 

513.0000 

-20,0000 

26,6000 

20.00 

6.2500 

.5000 

3.0300 

1.0100 

212.0000 

-21.1840 

26,6000 

06  ■ 

-4,500 
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OISPL.IWTS.)  » 

19875.00 

LONG,  C.8.  «  2,966 

<rwD. 

OF 

MIOSHIPS) 

OISPL.IVOL.)  ■ 

19862.03 

LONG,  C.G.  ■  3.212 

<r*D. 

OF 

MIOSHIPS) 

LONG,  GYRA0IU5  * 

116.489 

GM  »  3.722 
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Table  3.  Basic  Data  for  WOLVERINE  STATE,  Light  Load 


WOLVEHTME 

STATF.  LIGHT  LOAD. 

PAVIOSON 

LAB.  TEST 

CONDITION  - 

OCEANICS 

PROJ.  1093 

OPTION  CONTROL  tags 

•  AH 

C  0  F 

F  G  H 

1  J 

A  ? 

1  0  1 

0  0  1 

0  1 

NO. 

OF  STATIONS 

BASIC  INPUT  OATa 

length  * 

49ft#nn 

DENSITY  « 

028570 

OISPL.  « 

1210«5.on 

GRAVITY  *  32. 

175000 

STATION 

RtAM  A  OF  A  COEF. 

DRAFT 

Z-B'P 

HEIGHT 

ZFTa 

GYR.BOLL 

.SO 

9.750" 

.8*00 

17.3200 

7.7363 

551.0000 

-24.0610 

26. *500 

l.so 

34. soon 

.8700 

17.7600 

8.110* 

823,0000 

-10.0010 

26. *500 

2.50 

40.1000 

.9100 

18,2000 

8.55*0 

1670,0000 

-3.9210 

26.4500 

3.50 

70. 4000 

.9550 

18,6400 

9.0404 

1970,0000 

7,9990 

26.4500 

*.50 

71.5000 

.9750 

10.0800 

9.3810 

1067. OOLO 

11.4990 

24.4500 

5.50 

71.500" 

.9700 

19.5200 

9.5648 

1600,0000 

11.4990 

26.4500 

6.50 

71 .OOOn 

.9000 

19.9600 

9.31*7 

1390,0000 

7,4990 

26,4500 

7. SO 

62.A50O 

.7700 

2". *000 

8.4360 

1256,0000 

-1.3010 

28.4500 

8. SO 

37.000" 

.4900 

20.HAOO 

8.2665 

1077,0000 

-10,0010 

26.4500 

9.50 

8.400" 

.7*00 

21.2800 

8.7957 

685.0000 

-24.0010 

26.4500 

OG  *  3.200  GvRAn JUS .ROLL  *  26. *50 

CALCULATE  MOMENTS  at  STATION  5 


OERIVtn  RESULTS 

OISPL.IWTS.)  ■  1208Q.00 

LONG.  C.R.  «  .Tm  |FwD.  uF  MIDSHIPS)  OISPL.tVOL.)  •  121*7.21 

LONG.  C.G.  *  .401  IFwO.  OF  MIOSHIPS)  LONG.  GySAOIUS  »  123.391  GM  •  7,332 

Table  4.  Basic  Data  for  SERIES  60,  BLOCK  .80  Hull 

SEPIES  40  Mull  Foom,  0.80  Pluck  (TNO  RPT.  NO.  100  SI  OCEANICS  PROJECT  NO.  1093 

OPTION  CONTROL  TAGS  -APCOEFGHIJ 

1  2  l  3  1  0  ?  0  I  1  NO.  OF  STATIONS  »  ?0 

basic  input  data 


length  *  193.00  DENSITY  »  1.025000 

OISPL.  *  *0124.4"  GRAVITY  a  9.806650 


STATION 

BEAM 

AREA  COEF. 

DRAFT 

2-BAR 

HEIGHT 

ZETA 

GyR.ROLL 

0.00 

n.ooon 

0.0000 

0.0000 

0,0000 

240.6000 

0.0000 

8,9602 

1.00 

1A.3900 

.8720 

11.0300 

5.0*4* 

*81.3000 

o.oono 

8.9602 

2.00 

22.850" 

.89*0 

11.0300 

5.1253 

1203.2000 

0,0000 

8,9602 

3.00 

26,580" 

.9290 

11.0100 

5.254" 

2406.3000 

0,0000 

8.9602 

4,00 

27.540" 

.9700 

11.0300 

5.40*7 

3850.1000 

o.oono 

8.9602 

5.00 

27.5700 

.9910 

11.0300 

5,4819 

*090.7000 

o.oono 

8.9602 

6.00 

27.570" 

.99*0 

11.0300 

5.4929 

4331.4000 

0.0000 

8.9602 

7.00 

27.570" 

.99*0 

11.0300 

5.4929 

*331.4000 

O.oono 

8.9602 

8.00 

27.5700 

.99*0 

11.0300 

5.4929 

3368.8000 

0.0000 

8.9602 

9.00 

27.5700 

.99*0 

11.0300 

5.4929 

1684.4000 

0.0000 

8.9602 

10.00 

27.5700 

.9940 

11.0300 

5.4929 

1684.4000 

o.oono 

8.9602 

11.00 

27,5700 

.99*0 

11,0300 

5,4929 

1**3.8000 

o.oono 

8.9602 

12.00 

27.570" 

.9930 

11.0300 

5.4893 

2195,8000 

o.oono 

8.9602 

13.00 

27.570" 

.9890 

11,0300 

5.*7*4 

3290.7000 

0.0000 

8.9602 

14.00 

27,5700 

.9680 

11.0300 

5.3973 

3633.6000 

o.oono 

8.9602 

15.00 

27.240" 

.9210 

11.0300 

5.22*5 

3*65.1000 

0.0000 

8,9602 

16.00 

25.9*0" 

.8510 

11.0300 

4.9672 

31*6.3000 

0,0000 

8.9602 

17.00 

23.460" 

.7580 

11.0300 

4.4252 

1955.1000 

o.oono 

8.9602 

18.00 

19.630" 

.4270 

11.0300 

4.1436 

721.9000 

0,0000 

8.9602 

19.00 

13.870" 

.4190 

11.0300 

3.3789 

481.3000 

o.oono 

8.9602 

20.00 

4,4100 

.5300 

1.1000 

.3777 

120.3000 

o.oono 

8.9602 

OG  «  -1.099  GvPAntUS.HOLL  «  8.960 
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OERIVtO  results 


OISPL.IWTS.)  « 

*8124,50 

LONG.  C.G,  « 

*,714  (FhO. 

OF 

uIOSHlPSI 

DISPL.IVOL.)  ■ 

*8077.53 

LONG,  C.G,  * 

4. BPS  (FHO, 

OF 

MIOSHIPS) 

LONG.  GYRADIUS  ■ 

*6.159 

GM  ■ 

27 


Table  5.  Basic  Data  for  SERIES  60,  BLOCK  .70  Hull 


SEHIFS  *.0  HULL  FnPw,  .70  BLOCK,  L/fl  »  7.0,  L/H  •  17.5  (NSMR  TESTS)  OCEANICS  1093 

I 

wo.  or  stations  «  ?o 


OPTION  COWTHOL  TAOS  -  A 
1 


t>  c 

2  1 


E  F  0  M 
10  3  1 


BASIC  INPUT  DATA 


LENOTH 

OISPL. 


400,00 

Iiusa.oo 


DFN5ITY 

GRAVITY 


.038571 

33.175000 


STATION 

OtAM 

Aota  CoEF 

.  DRAFT 

z-bar 

WEIGHT 

ZFTA 

OVP.ROLL 

0.00 

.  1  700 

.QO00 

1 .0000 

.4667 

52.3000 

-9.10TS 

21.3500 

1.00 

9.4300 

.8486 

22.8600 

10,2763 

104.6000 

-6,6075 

21.3500 

2.00 

21 .»800 

.8341 

22.8600 

10.1658 

209.1000 

-4,6075 

21.3500 

3.00 

34. MOO 

.8619 

22.8600 

10.3777 

366.0000 

-2. 6075 

21.3500 

4.00 

44.9700 

.8908 

22.8600 

10.5979 

679,6000 

•1 . 1075 

21.3500 

5.00 

52, noon 

,9210 

22.8600 

10.8280 

815.6000 

-.0075 

21.3500 

4.00 

00,4000 

.9526 

22.8600 

11.0688 

1045.6000 

.8975 

21.3500 

7.00 

,9100 

.0761 

22.8600 

11.2479 

993,3000 

1.3905 

21.3500 

8.00 

57,1400 

.9850 

22.8600 

11.3157 

313.7000 

i.bops 

21.3500 

9.00 

57.1400 

,9860 

22.8600 

11.1233 

313.7000 

1,8975 

21.3500 

10,00 

57,1400 

,Q88U 

22,8600 

11.3233 

334.6000 

1.8925 

21.3500 

ll.no 

57.140« 

.9860 

22.0600 

11.1233 

444.4000 

1.8975 

21.3500 

12.00 

87.1400 

.0830 

22.6600 

11.3005 

627.4000 

1,8975 

21.3500 

13.  no 

57.0800 

.0672 

22.8600 

11.1801 

836.5000 

1.3975 

21.3500 

14.00 

54.8800 

.9363 

22.8600 

10.9446 

836.5000 

,8975 

21.3500 

15.00 

«8#^70ft 

.8812 

22,8600 

10.5247 

862.6000 

-,00T5 

21.3500 

l*.OQ 

52.1100 

•  8109 

22.8600 

9.9091 

784.2000 

-I.lOTS 

21.3500 

17.00 

48,8800 

.7133 

22.8600 

9.2453 

470.5000 

-2. 6075 

21.3500 

l«.nO 

18,4000 

.5869 

22.8600 

0,2822 

209.1000 

-*. 6075 

21.3500 

19.00 

28,4000 

,  1884 

22.8600 

6.7694 

104.6QC0 

-6.6075 

21.3500 

20,00 

8.0100 

.5060 

2.2H00 

.7646 

52.3000 

-9.1075 

21.3500 

00  • 

2.554 

GvRAOIUSb 

HULL  •  21 

.350 

CALCULATE  MOMENTS  AT  STATION  10 

OEPtvtn  ofkuLTS 


LONG.  C.“.  «  -7.**.-)  (FwD,  uF  HIOSHIPS) 
lcno.  r.o.  •  -s.no n  tfwn,  of  wlnSHlPS) 


OISPL.  (MTS.)  *  10*5N.?0 

OISPL. (VOL.)  »  10**4.73 

LONS.  STPAOIUS  »  95.080 


SM  ■ 


3.  *57 


Table  6.  Basic  Data  for  T-2  Tanker  Model 


T-S  TAN.FP  HOOFi  .  OFVIOSON  LAP.  (NUMATA)  TFST  CONDITIONS 


OPTION  COWTHOL  T»G« 
BASIC  INPUT  0*T* 


BCD 
3  0  0 


E  F  6  H 
10  3  1 


I  J 
0  1 


OCEANICS  PROJECT  1093 

NO.  OF  STATIONS  ■  30 


length  9 

4,00 

DENSITY 

»  62. 

500000 

OISPL.  * 

41  .0? 

G9AVITY 

«  32. 

1 75000 

STaTJOm 

PLAM 

AOEA  COEF.  DRAFT 

Z-8AP 

HEIGHT 

ZETa 

GYP. POLL 

ft. 00 

0,0000 

0.0000 

•  2860 

.0*77 

.1313 

o.oono 

.2550 

1.00 

.1670 

.8375 

•  2060 

.1275 

.32*3 

0,0000 

.2550 

?.oo 

.3450 

.0716 

,2060 

.1300 

1 .2278 

o.oooo 

.2550 

3.00 

.4870 

,0616 

•  2860 

.1295 

1.5598 

0.0000 

.2550 

4.00 

.5810 

.0967 

•  2860 

.1332 

2.32*3 

o.oooo 

.2550 

5.00 

♦  6330 

.9390 

.2860 

.1372 

2.59** 

0,0000 

.2550 

6.  AO 

.6400 

.9766 

,2860 

.1*08 

2.6795 

0,0000 

.2550 

7. no 

•  64M 

.9020 

•  2660 

.1*17 

2.872* 

0,0000 

.2550 

0.00 

,6400 

.9020 

•  2860 

.1*13 

3.1351 

0,0000 

.2550 

9. no 

,6480 

.9820 

,2860 

.1*13 

2.8108 

0,0000 

.2550 

lo.oo 

.6480 

.9020 

.2860 

.1*13 

2.7799 

0,0000 

.2550 

11. no 

•  640^ 

.9820 

•  2860 

.1*13 

2.7799 

0,0000 

.2550 

1J».0C 

,6480 

.9820 

•  2860 

.1*13 

2.7876 

0,0000 

.2550 

I3.ft0 

.6480 

.9659 

.2860 

.1397 

2.7568 

0.0000 

.2550 

14.00 

.6400 

.9309 

•  2060 

.137? 

2.7027 

o.oooo 

.2550 

16.00 

.6410 

.9000 

•  2860 

.1335 

2.2780 

0,0'  : 

.2550 

16.00 

,6140 

.0428 

•  2860 

.1280 

2.0386 

0.00'  7 

.2550 

17.0y 

•  5500 

.7438 

•  2860 

.118* 

1,37*5 

0,0000 

.2550 

18.00 

,4350 

.5707 

.2060 

.1028 

1.0*25 

o.oooo 

.2550 

19.00 

*2030 

•  4324 

•  2860 

.0889 

,*788 

o.oooo 

.2550 

20.00 

.1020 

.5100 

.0520 

.0175 

.3398 

0,0000 

.2550 

1 

X 

•j 

o 

.OTO 

GvPAOIUS.ROLL  « 

.255 

CALCULATE  MOMENTS  AT  STATION  10 

OEPTvtn  pf«ULTS 


OISPL. (NTS.)  * 


*1.03 


LONG.  f.B. 
LONG.  c.O. 


.03n  (FHO.  OF  MIOSMIPS) 
.03)  (F.O.  OF  MIOSMIPS) 


OISPL. (VOL.)  * 
LONS.  GYRA01US  • 


»n.9* 

1,090  GM 


.053 


HAy 


The  results  of  the  calculations  are  presented  in  the  same 
form  as  the  experimental  data  in  the  various  sources.  For  the 
Wolverine  State,  the  results  are  given  for  the  full-scale  ship. 
For  the  T-2  Tanker,  model-scale  results  are  shown.  For  the 
Series  60  hulls,  results  are  shown  in  non-dimensional  form,  as 
follows: 


1 


I 


i 


Froude  Number:  F  = 

n 

Non-dimensional  wave  frequency  = 


Non-dimensional  moment 


Non-dimensional  shear 


where  B*  =  waterline  beam  amidships. 

The  comparison  between  calculations  of  vertical  and  lateral 
bending  moments  and  the  experimental  results  for  the  Wolverine 
State,  shown  in  Figures  2  and  3,  indicates  generally  very  good 
agreement.  This  holds  for  both  loading  conditions,  both  speeds, 
and  over  the  range  of  wave  angle  and  wavelength.  The  experimental 
results  shown  for  lateral  bending  moment  in  head  and  following  seas, 
where  lateral  motions  and  loads  should  be  zero  as  in  the  calculations, 
are  regarded  as  indicative  of  the  possible  error,  or  range  of  dis¬ 
crepancy,  to  be  expected  between  calculations  and  experimental  re¬ 
sults.  These  loads  are  believed  to  arise  in  the  model  tests  due 
to  its  free-running,  but  rudder  controlled,  condition.  That  is, 
the  model  may  undergo  small  lateral  motions,  with  rudder  corrections 
to  keep  course,  which  leads  to  the  measured  lateral  bending  moments. 

The  comparison  for  the  Series  60,  block  0.80  hull  shown  in 
Figure  4  for  vertical  and  lateral  bending  moments  indicates  excellent 
agreement,  it'  general.  Similar  results  were  also  shown  for  this 
hull  by  Faltinsen  [21]  based  on  a  new  strip  theory  of  Salvesen,  et 
al  [22] .  Figure  5  shows  the  torsional  moment  comparison,  while 
in  Figure  6  the  vertical  and  lateral  shear  forces,  which  were  also 
measured  by  Wahab  [18] ,  are  shown.  The  agreement  for  torsional 
moments  is  only  fair  and  indicates  excessive  response  at  roll  re¬ 
sonance  conditions.  The  agreement  for  the  shear  forces  is  quite 
good,  in  general,  with  the  exception  of  some  deviation  in  lateral 
shear  at  110°  wave  angle.  However,  the  shear  forces  are  generally 
small  at  midships,  and  should  really  be  investigated  at  the  quarter- 
length  points.  Vertical  and  lateral  bending  moment  responses  in 
irregular  seas  are  shown  in  Figure  7.  The  experimental  results 
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Figures  2—12  are  grouped  at  end  of  raport  beginning  with  page  33. 
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are  not  from  direct  irregular  wave  testing,  but  rather  are  calcu¬ 
lated  from  the  regular  wave  unidirectional  data,  using  the  particular 
sea  spectrum  indicated.  The  difference  between  long  crested  and 
short  crested  seas  results  are  particularly  interesting  for  the 
lateral  bending  moment.  They  show  that  while  the  response  is  min¬ 
imal  in  unidirectional  beam  seas,  compared  to  the  peaks  at  bow 
and  stern  quartering  headings,  the  short  crested  seas  response  is 
maximum  in  beam  seas. 

Figures  8  and  11  show  the  comparison  for  vertical  and  lateral 
bending  moments  for  the  Series  60,  block  0.70  hull  form.  A  wide 
range  of  ship  speed  is  covered  in  this  data.  The  T-2  Tanker  model 
comparison  is  shown  in  Figure  12.  For  the  150  and  120  degree  wave 
angles,  experimental  data  and  calculations  are  shown  over  a  range 
of  3peed  for  two  wavelengths,  i.e.  a  wavelength  equal  to  model 
length  and  a  wavelength  such  that  its  "effective  length"  is  equal 
to  model  length.  In  the  latter  case,  the  actual  wavelength  equalled 
the  model  length  times  the  cosine  of  the  wave  angle.  This  data 
covers  vertical  and  lateral  bending,  and  torsional  moments.  In 
general,  the  agreement  is  fairly  satisfactory,  considering  the  factors 
involved  in  the  experimental  comparison.  With  regard  to  this  point, 
consider  the  double  peak  calculated  vertical  bending  moment  response 
for  the  T-2  Tanker  at  120°  wave  heading  and  1.65  fps  model  speed 
(Figure  12h) .  While  the  corresponding  experimental  data  do  not  in¬ 
dicate  such  a  response  similar  double  peaked  responses  for  vertical 
bending  are  confirmed  by  experimental  results  for  Wolverine  State, 
full  load  (Figure  2c),  and  the  Series  60,  block  0.80  hull  (Figure  4b). 
The  greater  resolution  of  the  test  data  due  to  testing  at  more  wave¬ 
length  conditions  for  these  latter  cases  tends  to  produce  such  re¬ 
sults,  thereby  limiting  the  utility  of  the  experimental  points  for 
the  T-2  Tanker  as  a  complete  measure  of  bending  moment  variation. 

The  preceding  comparisons  have  demonstrated  the  capability 
of  the  present  analysis  and  its  computer  implementation  to  provide 
valid  predictions  of  wave-induced  structural  loads  on  conventional 
ship  hull  forms.  As  discussed  previously,  the  technique  used  is 
based  upon  a  sectional  representation  with  Lewis  forms,  and  hence 
bulbous  bows  cannot  be  represented  accurately  (i.e.  in  matching  the 
desired  sectional  form  with  the  resulting  shape  obtained  by  the 
Lewis  form  fit) .  However  some  limited  results  obtained  by  comparing 
the  outputs  from  a  Lewis  form  representation  with  that  from  an 
accurate  "close  fit'  technique  (see  [23])  showed  little  effect  on 
the  resulting  motions  of  heave  and  pitch  when  using  either  method 
of  determining  the  two-dimensional  sectional  added  mass  and  damping, 
although  the  inability  to  match  the  section  form  was  demonstrated. 

This  result  would  appear  to  imply  that  the  use  of  the  Lewis  form 
fit  produced  sufficiently  useful  data  for  sectional  forces  that 
would  manifest  whatever  influence  was  exhibited  by  the  bulbous  bow 
form,  or  possibly  that  such  a  localized  force  did  not  have  a  signi¬ 
ficant  influence  on  the  overall  body  motions.  In  either  case  the 
same  characteristics  would  be  expected  to  carry  over  as  well  to  the 
case  of  the  computation  of  bending  moments,  and  hence  the  presently 
developed  technique  can  also  be  useu  for  predictions  for  the  case 
of  bulbous  bow  hulls.  Since  the  computation  of  the  sectional  added 
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mass  and  damping  are  determined  by  a  specific  subroutine  in  the 
overall  computer  program  in  [3] ,  and  only  a  limited  portion  of  the 
hull  (at  the  bow)  is  affected  by  the  bulb,  the  use  of  a  specialized 
procedure  for  that  region  can  be  adapted  if  desired,  based  on  the 
methods  and  computer  program  used  in  (9] ,  for  example,  or  any  other 
simple  computer  program  developed  to  encompass  bulbous  bow  hulls . 


CONCLUSIONS 

An  analytical  method  for  the  determination  cf  wave-induced 
moments  on  ships  has  been  developed,  implemented  (via  computer 
program) ,  and  successfully  evaluated  by  comparison  with  a  large 
body  of  model-scale  experimental  data.  It  should  prove  to  be  a 
valuable  aid  to,  and  integral  element  of,  the  fundamental  and 
rational  ship  structural  design  approach.  It  can  be  used  to 
predict  the  ship  motions  and  wave-induced  vertical  and  lateral 
bending  moments,  and  torsional  moment,  at  any  station  along  the 
length,  for  a  ship  traveling  at  any  heading  relative  to  long  or 
short  crested  seas. 

The  computer  program,  which  embodies  the  developed  method, 
is  documented  in  complete  detail  in  [33 .  It  can  be  used  in  the 
basic  ship  design  process  for  the  prediction  or  determination  of 
both  ship  motions  and  the  wave-induced  structural  loads.  The 
approach  and  implementation  are  straight-forward,  and  the  program 
is  efficient  in  regard  to  computer  time  usage. 

While  the  possibilities  for  use  of  the  analytical  method 
appear  great,  some  additional  development  work  would  seem  to  be 
in  order.  The  influence  of  rudder  effects  should  be  investigated. 
The  effect  of  the  rudder  and  control  system  upon  ship  motions  and 
loads  needs  some  careful  attention.  In  addition,  the  effects  due 
to  mean  drift  forces  and  moments,  manifested  by  leeway  angles  and 
mean  rudder  angles,  ought  to  be  determined.  The  present  evaluation 
of  the  method  indicates  that  such  effects  are  relatively  small, 
since  the  responses  of  interest  are  those  of  oscillatory  nature 
with  a  frequency  equal  to  the  encounter  frequency  in  regular  waves, 
but  a  fuller  understanding  of  their  influence  is  nevertheless  re¬ 
quired. 
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